Modeling DNA Structure, Elasticity and Deformations at the Base-pair Level 
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We present a generic model for DNA at tiie base-pair level. We use a variant of the Gay-Berne 
potential to represent the stacking energy between neighboring base-pairs. The sugar-phosphate 
backbones are taken into account by semi-rigid harmonic springs with a non-zero spring length. 
The competition of these two interactions and the introduction of a simple geometrical constraint 
leads to a stacked right-handed B-DNA-like conformation. The mapping of the presented model to 
the Marko-Siggia and the Stack-of-Plates model enables us to optimize the free model parameters so 
as to reproduce the experimentally known observables such as persistence lengths, mean and mean 
squared base-pair step parameters. For the optimized model parameters we measured the critical 
force where the transition from B- to S-DNA occurs to be approximately 140pN. We observe an 
overstretched S-DNA conformation with highly inclined bases that partially preserves the stacking 
of successive base-pairs. 
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PACS numbers: 87.14.Gg,87.15.Aa,87.15.La,61.41.+e 



I. INTRODUCTION 

FoUowingthe discovery of the double helix by Watson 
and Crick the structure and elasticity of DNA has 
been investigated on various length scales. X-ray diffrac- 
tion studies of single crystals of DNA oligomers have led 
to a detailed picture of possible DNA conformations 0,01 
with atomistic resolution. Information on the behavior 
of DNA on larger scales is accessible through NMR 
and various optical methods 0, 0| , video and electron 
microscopy An interesting development of the last 
decade are nanomechanical experiments with individual 
DNA molecules H E El El El which, for example, re- 
veal the intricate interplay of supercoiling on large length 
scales and local denaturation of the double-helical struc- 
ture. 

Experimental results are usually rationalized in the 
framework of two types of models: base-pair steps and 
variants of the continuum elastic worm-like chain. The 
first, more local, approach describes the relative loca- 
tion and orientation of neighboring base pairs in terms 
of intuitive par ameters such as twist, rise, slide, roll 
etc. 0, ITsl Il6l Ir^ . In particular, it provides a me- 
chanical interpretation of the biological function of par- 
ticular sequences E3- The second approach models 
DNA on length scales beyond the helical pitch as a 
worm-like chain (WLC) with empirical parameters de- 
scribing the resistance to bending, twisting and stretch- 
ing |20| . The results are in remarkable agree- 
ment with the nanomechanical experiments mentioned 
above 0| . WLC models are commonly used in order to 
address biologically important phenomena such as super- 
coilin g [23. 123. or the wrapping of DNA around his- 
tones [23. In principle, the two descriptions of DNA are 
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linked by a systematic coarse-graining procedure: From 
given (average) values of rise, twist, slide etc. one can 
reconstruct the shape of the corresponding helix on large 
scales ^3 |26l . Similarly, the elastic constant char- 
acterizing the continuum model are related to the local 
elastic energies in a stack-of-plates model p^ . 

Difficulties are encountered in situations which can- 
not be described by a linear response analysis around 
the undisturbed (B-DNA) ground state. This situation 
arises routinely during cellular processes and is there- 
fore of considerable biological interest jl^ . A character- 
istic feature, observed in many nanomechanical experi- 
ments, is the occurrence of plateaus in force-elongation 
curves 0, O, 0. These plateaus are interpreted as 
structural transitions between microscopically distinct 
states. While atomistic simulations have played an im- 
portant role in identifying possible local structures such 
as S- and P-DNA jll E3i t;his approach is hmited to 
relatively short DNA segments containing several dozen 
base pairs. The behavior of longer chains is interpreted 
on the basis of stack-of-plates models with step-type de- 
pendent parameters and free energy penalties for non-B 
steps. Realistic force-elongation are obtained by a suit- 
able choice of parameters and as the consequence of con- 
straints for the total extension and twist (or their con- 
jugate forces) I2I. Similar models describing the non- 
linear response of B-DNA to stretching or untwist- 
ing [30I Isij l predict stability thresholds for B-DNA due 
to a combination of more realistic, short-range interac- 
tion potentials for rise with twist-rise coupling enforced 
by the sugar-phosphate backbones. 

Clearly, the agreement with experimental data will in- 
crease with the amount of details which is faithfully rep- 
resented in a DNA model. However, there is strong ev- 
idence both from atomistic simulations [s^ as well as 
from the analysis of oligomer crystal structures that 
the base-pair level provides a sensible compromise be- 
tween conceptual simplicity, computational cost and de- 
gree of reality. While Lavery et al. j32|| have shown that 



Typeset by REVTeX 



2 



the base-pairs effectively behave a s rig id entities, the re- 
sults of El Hassan and Calladine [33 and of Hunter et 
al. 0, 1^3 suggest that the dinucleotide parameters ob- 
served in oligomer crystals can be understood as a con- 
sequence of van-der-Waals and electrostatic interactions 
between the neighboring base-pairs and constraints im- 
posed by the sugar-phosphate backbone. 

The purpose of the present paper is the introduction of 
a class of "DNA- like" -molecules with simplified interac- 
tions resolved at the base or base pair level. In order to 
represent the stacking interactions between neighboring 
bases (base pairs) we use a variant .36] of the Gay-Berne 
potential ^3 used in studies of discotic liquid crystals. 
The sugar-phosphate backbones are reduced to semi-rigid 
springs connecting the edges of the disks/ellipsoids. Us- 
ing Monte-Carlo simulations we explore the local stack- 
ing and the global helical properties as a function of 
the model parameters. In particular, we measure the 
effective parameters needed to describe our systems in 
terms of stack-of-plates (SOP) and worm-like chain mod- 
els respectively. This allows us to construct DNA models 
which faithfully represent the equilibrium structure, fluc- 
tuations and linear response. At the same time we pre- 
serve the possibility of local structural transitions, e.g. 
in response to external forces. 

The paper is organized as follows. In the second section 
we introduce the base-pair parameters to discuss the helix 
geometry in terms of these variables. Furthermore we dis- 
cuss how to translate the base-pair parameters in macro- 
scopic variables such as bending and torsional rigidity. 
In the third section we introduce the model and discuss 
the methods (MC simulation, energy minimization) that 
we use to explore its behavior. In the fourth section we 
present the resulting equilibrium structures, the persis- 
tence lengths as a function of the model parameters, and 
the behavior under stretching. 



II. THEORETICAL BACKGROUND 
A. Helix geometry 

To resolve and interpret X-ray diffraction studies on 
DNA oligomers the relative position and orientation of 
successive base-pairs are analyzed in terms of Rise (Ri) , 
Slide (SI), Shift (Sh), Twist (Tw), Roll (Ro), and Tih 
(Ti) jSg (see Fig. In order to illustrate the relation 
between these local parameters and the overall shape of 
the resulting helix we discuss a simple geometrical model 
where DNA is viewed as a twisted ladder where all bars 
lie in one plane. For vanishing bending angles with Ro = 
Ti = each step is characterized by four parameters: 
Ri, SI, Sh, and Tw 18]. Within the given geometry a 
base pair can be characterized by its position r and the 
angle of its main axis with the n/b-axis (n points into the 
direction of the large axis, b points into the direction of 
the small axis, and t, representing the tangent vector of 
the resulting helix, is perpendicular to the n-b- plane as it 
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FIG. 1: (Color online) Illustration of all six base-pair param- 
eters and the corresponding coordinate system. 



is illustrated in Fig. ^ . At each step the center points are 

displaced by a distance vS? + St? in the n — b— plane. 
The angle between successive steps is equal to the twist 
angle and the center points are located on a helix with 

radius r = ysPTsi?/(2 sin(Tw/2)). 

In the following we study the consequences of impos- 
ing a simple constraint on the bond lengths li and I2 rep- 
resenting the two sugar phosphate backbones (the rigid 
bonds connect the right and left edges of the bars along 
the n-axis respectively). Ri is the typical height of a 
step which we will try to impose on the grounds that it 
represents the preferred stacking distance of neighboring 
base pairs. We choose Ri = 3.3A corresponding to the 
B-DNA value. One possibility to fulfill the constraint 
li — I2 ~ I = 6 A is pure twist. In this case a relationship 
of the twist angle and the width of the base-pairs d, the 
backbone length I and the imposed rise is obtained: 



Tw = arccos 



2P + 2Ri^ 



(1) 



Another possibility is to keep the rotational orientation 
of the base pair (Tw = 0), but to displace its center in 
the n-b-plane, in which case Ri2 + Sl^ + Sh2 = /2. With 
Sh — 0, it results in a skewed ladder with skew angle 
arcsin(Sl/0/7r [3. 

The general case can be solved as well. In a first step 
a general condition is obtained that needs to be fulfilled 
by any combination of Sh, SI, and Tw independently of 
Ri. For non-vanishing Tw this yields a relation between 
Sh and SI: 



tan(Tw) = 



Sh 

sr 



(2) 



Using Eq. ((2Jl the general equation can finally be solved: 
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FIG. 2: (Color online) Illustration of DNA geometry for a 
diameter of d = 16A: (1) Twisted ladder with SI = Sh = 0, 
Ri = 3.3A, Tw ^ 27r/10, (2) Skewed ladder with Tw = Sh = 
0, Ri = 3.4A, SI « 5.0A, (3) Helix with Tw = 27r/12, Ri = 
3.4A, SlRi 2.7A, ShRi 1.6A. 



Sl = 



V2 



cos( — 



sec(^)2(2?2 



- Ri^) 



(3) 



Eq. Q is the result of the mechanical coupling of slide, 
shift and twist due to the backbones. Treating the rise 
again as a constraint the twist is reduced for increasing 
slide or shift motion. The center-center distance c of two 
neighboring base-pairs is given by 



c= \/Ri 



Sl^ (1 



tan(Tw)2). 



(4) 



For Tw = and a given value of Ri the center-center 
distance is equal to the backbone length I and for Tw = 
arccos ({(P - 2P + 2Ri^)/(f) one obtains c = Ri. 



B. Thermal fluctuations 

In this section we discuss how to calculate the effective 
coupling constants of a harmonic system valid within lin- 
ear response theory describing the couplings of the base- 
pair parameters along the chain. Furthermore we show 
how to translate measured mean and mean squared val- 
ues of the 6 microscopic base-pair parameters into macro- 
scopic observables such as bending and torsional persis- 
tence length. This provides the linkage between the two 
descriptions: WLC (worm-like chain) versus SOP (stack- 
of-plates) model. 

Within linear response theory it should be possible to 
map our model onto a Gaussian system where all trans- 
lational and rotational degrees of freedom are harmon- 
ically coupled. We refer to this model as the stack-of- 
plates (SOP) model The effective coupling constants 
are given by the second derivatives of the free energy in 



terms of base-pair variables around the equilibrium con- 
figuration. This yields 6x6 matrices /C"™ describing 
the couplings of the base-pair parameters of neighboring 
base-pairs along the chain: 



(5) 



Therefore one can calculate the (N — 1) x (iV — 1) cor- 
relation matrix C in terms of base-pair parameters. N is 
thereby the number of base-pairs. 



(C) = 



/C" /Ci2 K}^ K}^ 

_ I K?^ 



(6) 



The inversion of C results in a generalized connectivity 
matrix with effective coupling constants as entries. 

The following considerations are based on the assump- 
tion that one only deals with nearest-neighbor interac- 
tions. Then successive base-pair steps are independent 
of each other and the calculation of the orientational 
correlation matrix becomes feasible. In the absence of 
spontaneous displacements (SI = Sh = 0) and sponta- 
neous bending angles (Ti = Ro = 0) as it is the case for 
B-DNA going from one base-pair to the neighboring im- 
plies three operations. In order to be independent of the 
reference base pair one first rotates the respective base 
pair into the mid-frame with TZ{T'Wsp/i) {TZ is a rotation 
matrix, Tw^p denotes the spontaneous twist), followed 
by a subsequent overall rotation in the mid-frame 



A 




ti • bi+i ti • rii+i 
hi • bi+i bi • ni+i 
rii • bi+i rii • rii+i, 



(7) 



taken into account the thermal motion of Ro, Ti and 
Tw, and a final rotation due to the spontaneous twist 
TZ{T'Wsp/2). The orientational correlation matrix be- 
tween two neighboring base pairs can be written as 
{0,i+i) = n{Twsp/2){A)n{Twsp/2). A describes the 
fiuctuations around the mean values. As a consequence 
of the independence of successive base-pair parameters 
one finds (0,j) = iJliTwsp/2) {A) n.{Twsp/2)y~' where 
the matrix product is carried out in the eigenvector ba- 
sis of TZ{Twsp/2) {A) TZ[Twsp/2). In the end one finds a 
relationship of the mean and mean squared local base- 
pair parameters and the bending and torsional persis- 
tence length. The calculation yields an exponentially de- 
caying tangent-tangent correlation function (t(O)-t(s)) = 
exp(— s/^p) with a bending persistence length 



2(Ri) 



((Ti2) + (Ro^))" 



(8) 



In the following we will calculate the torsional per- 
sistence length. Making use of a simple relationship be- 
tween the local twist and the base-pair orientations turns 
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out to be more convenient than the transfer matrix ap- 
proach. 

The (bi)normal-(bi)normal correlation function is an 
exponentially decaying function with an oscillating term 
depending on the helical repeat length h — p(Ri) and 
the helical pitch p = 27r/(Tw) respectively, namely 
(n(0) • n(s)) = exp(— s/^„) cos(27r s//i). The torsional 
persistence length Z„ = 4 can be calculated in the fol- 
lowing way. It can be shown that the twist angle Tw 
of two successive base-pairs is related to the orientations 
{t,b, n} and {t',b',n'} through 



cos(Tw) 



n • n 



bb' 



1 + t • t' 



(9) 



Taking the mean and using the fact that the orientational 
correlation functions and twist correlation function decay 
exponentially 



exp{-l/lTn 



2exp(-l/?„) 
1 + exp(-l/Zp) 



(10) 



yields in the case of stiff filaments a simple expression of 
In depending on Ip and Itw- 
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~2 



h 
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where the twist persistence length is defined as 

(Ri) 



(T. 



(11) 



(12) 



III. MODEL AND METHODS 

Qualitatively the geometrical considerations suggest a 
B-DNA like ground state and the transition to a skewed 
ladder conformation under the influence of a sufficiently 
high stretching force, because this provides the possibil- 
ity to lengthen the chain and to partially conserve stack- 
ing. Quantitative modeling requires the specification of 
a Hamiltonian. 



A. Introduction of the Hamiltonian 

The observed conformation of a dinucleotide base-pair 
step represents a compromise between (i) the base stack- 
ing interactions (bases are hydrophobic and the base- 
pairs can exclude water by closing the gap in between 
them) and (ii) the preferred backbone conformation 
(the equilibrium backbone length restricts the conforma- 
tional space accessible to the base-pairs) [s^. Packer 
and Hunter have shown that roll, tilt and rise are 
backbone-independent parameters. They depend mainly 
on the stacking interaction of successive base-pairs. In 
contrast twist is solely controlled by the constraints im- 
posed by a rigid backbone. Slide and shift are sequence- 
dependent. While it is possible to introduce sequence 
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FIG. 3: (Color online) (left) Illustration of the underlying 
idea. The base-pairs are represented as rigid ellipsoids. The 
sugar-phosphate backbone is treated as semi-rigid springs 
connecting the edges of the ellipsoid, (right) Introduced in- 
teractions lead to a right-handed twisted structure. 



dependant effects into our model, they are ignored in the 
present paper. 

In the present paper we propose a generic model for 
DNA where the molecule is described as a stack of thin, 
rigid ellipsoids representing the base pairs (Fig. O. The 
shape of the ellipsoids is given by three radii a, b, c of 
the main axes in the body frames which can be used to 
define a structure matrix 



S 




(13) 



2a corresponds to the thickness, 2b to the depth which 
is a free parameter in the model, and 2c — ISA to the 
width of the ellipsoid which is fixed to the diameter of a 
B-DNA helix. The thickness 2a will be chosen in such a 
way that the minimum center-center distance for perfect 
stacking reproduces the experimentally known value of 
3.3A. 

The attraction and the excluded volume between the 
base pairs is modeled by a variant of the Gay-Berne po- 
tential [3^ 133 for ellipsoids of arbitrary shape Si , rela- 
tive position ri2 and orientation A^. The potential can 
be written as a product of three terms: 

t/(Ai, A2,fi2) = C/r(Ai, A2,fi2) 

X ?7i2(Ai, A2, fi2) Xl2(Ai, A2, fi2). 

(14) 

The first term controls the distance dependence of the 
interaction and has the form of a simple LJ potential 



Ur = 4eGB 



12 



h + 7(7 



(15) 



where the interparticle distance r is replaced by the dis- 
tance h of closest approach between the two bodies: 



min(|fi - fj\) 



(16) 



with i G Body 1 and j € Body 2. The range of interac- 
tion is controlled by an atomistic length scale a = 3. 3 A, 
representing the effective diameter of a base-pair. 
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In general, the calculation of h is non-trivial. We use 
the following approximative calculation scheme which is 
usually employed in connection with the Gay-Berne po- 
tential: 



/l(Ai, A2,ri2) = ri2 - cri2(Ai, A2,fi2) 



(17) 



ai2(Ai,A2,fi2) = [-fi^2Gr2'(Ai,A2)fi2]-^/^(18) 



Gi2(Ai,A2) = AfS^Ai + Afs2A2 



(19) 



In the present case of oblate objects with rather perfect 
stacking behavior Eq. produces only small deviations 
from the exact solution of Eq. (|16|l . 

The other two terms in Eq. (|14|l control the interaction 
strength as a function of the relative orientation A* A2 
and position ri2 of interacting ellipsoids: 



r7i2(Ai, A2,fi2) 



det[Si]/gj -|-det[S2]/cri 

(det[Hi2]/(ai -^(72)) 



1/2 



(20) 



Hi2(Ai,A2,fi2) - — AfS?Ai-f — Afs2A2(21) 

(Tl (72 



a^{A^,rl2) = {rl2 ^Ai fi2 



-1/2 



(22) 



and 



Xi2(Ai,A2,fi2) = [2ff2 Br2'(Ai,A2) fi2] (23) 
Bi2(Ai,A2) = AfEiAi+AjE2A2 (24) 



with 



E,; = a 



ai 
bi Ci 











bi 

di Ci 











o-i hi 



det[S, 



(25) 



We neglect electrostatic interactions between neighbor- 
ing base-pairs since at physiological conditions the stack- 
ing interaction dominates 0, • 

At this point we have to find appropriate values for 
the thickness 2a and the parameter 7 of Eq. H15|) . Both 
parameters influence the minimum of the Gay-Berne po- 
tential. There are essentially two possible procedures. 
One way is to make use of the parameterization result 
of Everaers and Ejtehadi Hil, i.e. 7 = 2^/^ - 30-l/^ 
and to choose a value of a ~ 0.7 that yields the mini- 
mum center-center distance of 3. 3 A for perfect stacking. 
Unfortunately it turns out that the fluctuations of the 
bending angles strongly depend on the flatness of the el- 
lipsoids. The more flat the ellipsoids are the smaller are 
the fluctuations of the bending angles so that one ends up 
with extremely stiff filaments with a persistence length 
of a few thousand base-pairs. This can be seen clearly 
for the extreme case of two perfectly stacked plates: each 
bending move leads then to an immediate overlap of the 
plates. That is why we choose the second possibility. We 
keep 7 as a free parameter that is used in the end to 
shift the potential minimum to the desired value and fix 
the width of the ellipsoids to be approximately half the 
known rise value a = 1.55A. This requires 7 = 1.07. 



The sugar phosphate backbone is known to be nearly 
inextensible. The distance between adjacent sugars 
varies from 5.5A to 6.5A 18] . This is taken into account 
by two stiff springs with length h — I2 — 6.0A connect- 
ing neighboring ellipsoids (see Fig. The anchor points 
are situated along the centerline in n-direction (compare 
Fig. Hand Fig. O with a distance of ±8A from the center 
of mass. The backbone is thus represented by an elastic 
spring with non-zero spring length ~ SA 



He 



k 

2 



|ri,j+i - ri,,,| - lof + (|r2,,+i - r2,j| - ^o)^] . 

(26) 

Certainly a situation where the backbones are brought 
closer to one side of the ellipsoid so as to create a mi- 
nor and major groove would be a better description of 
the B-DNA structure. But it turns out that due to the 
ellipsoidal shape of the base-pairs and due to the fact 
that the internal base-pair degrees of freedom (propeller 
twist, etc.) cannot relax a non-B-DNA-like ground state 
is obtained where roll and slide motion is involved. 

The competition between the GB potential that forces 
the ellipsoids to maximize the contact area and the har- 
monic springs with non-zero spring length that does not 
like to be compressed leads to a twist in either direction 
of the order of ±7r/5. The right-handedness of the DNA 
helix is due to excluded volume interactions between the 
bases and the backbone [l^ which we do not represent 
explicitly. Rather we break the symmetry by rejecting 
moves which lead to local twist smaller than — 7r/18. 

Thus we are left with three free parameters in our 
model, the GB energy depth e = min(J7) which con- 
trols the stacking interaction, the spring constant k which 
controls the torsional rigidity, and the depth b of the el- 
lipsoids which influences mainly the fluctuations of the 
bending angles. All other parameters such as the width 
and the height of the ellipsoids, or the range of interac- 
tion a = 3.3A which determines the width of the GB 
potential are fixed so as to reproduce the experimental 
values for B-DNA. 



B. MC simulation 

In our model all interactions are local and it can there- 
fore conveniently be studied using a MC scheme. In addi- 
tion to trial moves consisting of local displacements and 
rotations of one ellipsoid by a small amplitude, it is possi- 
ble to employ global moves which modify the position and 
the orientation of large parts of the chain. The moves are 
analogous of (i) the well-known pivot move , and (ii) a 
crankshaft move where two randomly chosen points along 
the chain define the axis of rotation around which the in- 
ner part of the chain is rotated. The moves are accepted 
or rejected according to the Metropolis scheme |4l] |. 

Fig. 21 shows that these global moves significantly im- 
prove the efficiency of the simulation. We measured the 
correlation time r of the scalar product of the tangent 
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FIG. 4: (Color online) Time correlation functions of the 
scalar product of the tangent vectors of the first and the 
last monomer r = ^(0,1) • t{t,N) with iV = 10 (red), A'' = 
20 (green), A'' = 50 (blue) for (a) global and (b) local moves. 
It is observed that Tgiobai is independent of the chain length 
whereas Tiocai scales as N^. The 'time' is measured in units 
of sweeps where one MC sweep corresponds to A'^ trials. The 
CPU time for one sweep scales as A^^ in case of global moves 
and as A'^ in case of local moves. Thus the simulation time t 
scales as tiocai oc A^** and tgiobai oc A^^. 



vectors of the first and the last monomer of 200 inde- 
pendent simulation runs with N = 10, 20, 50 monomers 
using (i) only local moves and (ii) local and global moves 
(ratio 1:1). The correlation time of the global moves is 
independent of the chain length with Tgiobai ~ 78 sweeps 
whereas Tiocai scales as N^. 

Each simulation run comprises 10^ MC sweeps where 
one MC sweep corresponds to 2N trials (one rotational 
and one translational move per base pair) with N denot- 
ing the number of monomers. The amplitude is chosen 
such that the acceptance rate equals approximately to 
50%. Every 1000 sweeps we store a snapshot of the DNA 
conformation. We measured the 'time' correlation func- 
tions of the end-to-end distance, the rise of one base- 
pair inside the chain and all three orientational angles 



of the first and the last monomer and of two neighbor- 
ing monomers inside the chain in order to extract the 
longest relaxation time T,nax- We observe Tmax < 1000 
for all simulation runs. 

An estimate for the CPU time required for one sweep 
for chains of length = 100 on a AMD Athlon MP 
2000-1- processor results in 0.026s which is equivalent to 
*s per move. 



1.33 lO"'^^ 



C. Energy minimization 

We complemented the simulation study by zero tem- 
perature considerations that help to discuss the geomet- 
ric structure that is obtained by the introduced interac- 
tions and to rationalize the MC simulation data. Fur- 
thermore they can be used to obtain an estimate of the 
critical force fcrit that must be applied to enable the 
structural transition from B-DNA to the overstretched 
S-DNA configuration as a function of the model param- 
eters {e, k, b}. 



IV. RESULTS 

In the following we will try to motivate an appropriate 
parameter set {e, fc, b} that can be used for further inves- 
tigations within the framework of the presented model. 
Therefore we explore the parameter dependence of ex- 
perimental observables such as the bending persistence 
length of B-DNA Ij, » 150bp, the torsional persistence 
length It « 260bp [42|, the mean values and correlations 
of all six base-pair parameters and the critical pulling 
force fcrit ~ 65pN jll|, i43|, (4^, |45i that must be apphed 
to enable the structural transition from B-DNA to the 
overstretched S-DNA configuration. In fact, static and 
dynamic contributions to the bending persistence length 
Ip of DNA are still under discussion. It is known that 
Ip depends on both the intrinsic curvature of the double 
helix due to spontaneous bending of particular base-pair 
sequences and the thermal fluctuations of the bending 
angles. Bensimon et al. [46] introduced disorder into 
the WLC model by an additional set of preferred ran- 
dom orientation between successive segments and found 
the following relationship between the pure persistence 



length Ip 



i.e. without disorder, the effective per- 



sistence length leff and the persistence length Idisorder 
caused by disorder: 



''pure 



1 - 



< 1 
> 1 



(27) 



Since we are dealing with intrinsically straight filaments 
with 1/ Idisorder = 0, wc mcasure Imire- Rcceut esti- 
mates of Idisorder range between 430 Wl\ and 4800 
base-pairs using cryo-electron microscopy and cyclization 
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(Ri) 


(Sh) 


(Sl> 


(Tw> 


(Ti) 


{Ro> 


(c) 


Ip 





3.26 


0.0 


0.0 


0.64 


0.0 


0.0 


3.26 


oo 


1 


3.37 


0.01 


-0.01 


0.62 


0.0 


0.0 


3.47 


172.8 


2 


3.76 


-0.01 


-0.03 


0.47 


0.0 


0.0 


4.41 


25.3 


3 


4.10 


-0.01 


0.01 


0.34 


0.0 


-0.01 


5.07 


14.4 


5 


4.30 


0.03 


-0.02 


0.27 


0.0 


0.01 


5.39 


13.6 



TABLE I: Dependence of mean values of all six step param- 
eters and of the mean center-center distance (c) on the tem- 
perature for 2b = llA, e = 20fcsr, k = GAksT/A^^ (Ri), 
(Sh), (SI) and (c) are measured in [A], Ip in base-pairs. 

experiments respectively implicating values between 105 
and 140 base-pairs for Ipure- 

A. Equilibrium structure 

As a first step we study the equilibrium structure of 
our chains as a function of the model parameters. To 
investigate the ground state conformation we rationalize 
the MC simulation results with the help of the geometri- 
cal considerations and minimum energy calculations. In 
the end we will choose parameters for which our model 
reproduces the experimental values of B-DNA • 

(Ri) = 3.3-3.4A 

(SI) = oA 

(Sh) = OA 

(Tw) = 27r/10.5- 27r/10 
(Ti) - 
(Ro) = 0. 

We use the following reduced units in our calculations. 
The energy is measured in units of ksT , lengths in units 
of A, forces in units of UbTIC^ « 40pN. 

We start by minimizing the energy for the various con- 
formations shown in Fig. |21 to verify that our model 
Hamiltonian indeed prefers the B-Form. Since we have 
only local (nearest neighbor) interactions we can restrict 
the calculations to two base pairs. There are three lo- 
cal minima which have to be considered: (i) a stacked, 
twisted conformation with Ri = 3.3, SI, Sh, Ti, Ro — 
0, Tw = 7r/10, (ii) a skewed ladder with Ri = 3.3, SI = 
5.0, Sh, Tw, Ti, Ro = 0, and (iii) an unwound helix with 
Ri = 6.0, SI, Sh, Ti, Ro = 0, Tw = 0. Without an exter- 
nal pulling force the global minimum is found to be the 
stacked twisted conformation. 

We investigated the dependence of Ri and Tw on the 
GB energy depth e that controls the stacking energy for 
different spring constants k. Ri depends neither on e nor 
on k nor on b. It shows a constant value of Ri « 3.3A for 
all parameter sets {e, fc, 6}. The resulting Tw of the mini- 
mum energy calculation coincides with the geometrically 
determined value under the assumption of fixed Ri up to 
a critical e. Up to that value the springs behave effec- 
tively as rigid rods. The critical e is determined by the 
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FIG. 5: (Color online) (a) Rise [A] and (b) twist as a func- 
tion of e [fcsT] for 26 = 8, 9, 10, llA (red, green, blue, 
purple). For every h there are two data sets for k — 
32 (plus), 64 (circles) [/cbT/A^]. The dotted line corresponds 
to the minimum energy value. (Ri) depends only on e. In the 
limit of e — > oo the minimum energy value is reached, 
(b) In addition to the MC data and the minimum energy cal- 
culation we calculated the twist with Eq. Q using the mea- 
sured mean rise values of (a). One can observe that (Tw) 
changes with all three model parameters. Increasing y and 
k decreases especially the fluctuations of Tw and Sh so that 
(Tw) increases as a result of the mechanical coupling of the 
shift and twist motion. In the limit of e, fc — > cxa the minimum 
energy value is reached. 

torque T{k,e) that has to be applied to open the twisted 
structure for a given value of Ri. 

Using MC simulations we can study the effects arising 
from thermal fluctuations. Plotting (Ri), and (Tw) as a 
function of the GB energy depth e one recognizes that 
in general (Ri) is larger than Ri(T = 0). It converges 
only for large values of e to the minimum energy values. 
This can be understood as follows. Without fluctuations 
the two base pairs are perfectly stacked taking the mini- 
mum energy configuration Ri = 3.3A, SI, Sh, Ti, Ro — 0, 
and Tw = tt/IO. As the temperature is increased the 
fluctuations can only occur to larger Ri values due to 
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the repulsion of neighboring base pairs. A decrease of Ri 
would cause the base-pairs to intersect. Increasing the 
stacking energy reduces the fluctuations in the direction 
of the tangent vector and leads to smaller (Ri) value. In 
the limit e — *■ oo it should reach the minimum energy 
value which is observed from the simulation data. In 
turn the increase of the mean value of rise results in a 
smaller twist angle (Tw) . We can calculate with the help 
of Eq. (Q) the expected twist using the measured mean 
values of (Ri). Fig. shows that there is no agreement. 
The deviations are due to fluctuations in SI and Sh which 
cause the base-pairs to untwist. This is the mechanical 
coupling of SI, Sh, and Tw due to the backbones already 
mentioned in section III Al It is observed that a stiffer 
spring k and a larger depth of the ellipsoids b result in 
larger mean twist values. Increasing the spring constant 
k means decreasing the fluctuations of the twist and, due 
to the mechanical coupling, of the shift motion around 
the mean values which explains the larger mean twist 
values. An increase of the ellipsoidal depth b in turn de- 
creases the fluctuations of the bending angles. The cou- 
pling of the tilt fluctuations with the shift fluctuations 
leads to larger values for (Tw). The corresponding limit 
where (Tw) — > Tw(T = 0) is given by fc, e ^ oo. 

The measurement of the mean values of all six base- 
pair step parameters for different temperatures is shown 
in Table One can see that with increasing tempera- 
ture the twist angles decrease while the mean value of rise 
increase. The increase of the center-center distance is not 
only due to fluctuations in Ri but also due to fluctuations 
in SI and Sh. That is why there are strong deviations of 
(c) from (Ri) even though the mean values of SI and Sh 
vanish. Note that the mean backbone length (l) always 
amounts to about 6A. 

The calculation of the probability distribution func- 
tions of all six base-pair parameters shows that especially 
the rise and twist motion do not follow a Gaussian behav- 
ior. The deviation of the distribution functions from the 
Gaussian shape depends mainly on the stacking energy 
determined by e. For smaller values of e one observes 
larger deviations than for large e values. 

It is worthwhile to mention that there are mainly two 
correlations between the base-pair parameters. The first 
is a microscopic twist-stretch coupling determined by a 
correlation of Ri and Tw, i.e. an untwisting of the he- 
lix implicates larger rise values. A twist-stretch coupling 
was introduced in earlier rod models ji^ [sO, E| moti- 
vated by experiments with torsionally constrained DNA 
|52| which allow for the determination of this constant. 
Here it is the result of the preferred stacking of neighbor- 
ing base-pairs and the rigid backbones. The second cor- 
relation is due to constrained tilt motion. If we return to 
our geometrical ladder model we recognize immediately 
that a tilt motion alone will always violate the constraint 
of fixed backbone length I. Even though we allow for 
backbone fiuctuations in the simulation the bonds are 
very rigid which makes tilting energetically unfavorable. 
To circumvent this constraint tilting always involves a 
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FIG. 6; (Color online) Contour plots of measured clouds for 
rise- twist, shift-tilt, and roll-tilt to demonstrate internal cou- 
plings and the anisotropy of the bending angles (26 = llA, 
e = 20fcsr, k = 64kBT/A^). 



directed shift motion. 

Fig. shows that we recover the anisotropy of the 
bending angles Ro and Ti as a result of the spatial di- 
mensions of the ellipsoids. Since the overlap of successive 
ellipsoids is larger in case of rolling it is more favorable 
to roll than to tilt. 

The correlations can be quantified by calculating the 
correlation matrix C of Eq. Inverting C yields the 
effective coupling constants of the SOP model K. = C~^. 
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FIG. 7: (Color online) Comparison of probability distribu- 
tion functions of all base-pair parameters for e — 20kBT, 
k = 64fcsT/A^, 26 = 8 A. The Gaussians are plotted with 
the measured mean and mean squared values of the MC sim- 
ulation. 



Due to the local interactions it suffices to calculate mean 
and mean squared values of Ri, SI, Sh, Tw, Ro, and Ti 
characterizing the 'internal' couplings of the base-pairs: 

C = (a),,, V*,j €{!,..., 6} (28) 

with a^cy = (xy) - {x)(y). 



B. Bending and torsional rigidity 

The correlation matrix of Eq. (|28|l can also be used to 
check eqs. © and Hll|) . Therefore we measured the ori- 
entational correlation functions (ti-tj), {rii-nj), (hi-hj) 
and compared the results to the analytical expressions as 
it is illustrated in Fig. |S1 The agreement is excellent. 

The simulation data show that the bending persistence 
length does not depend on the spring constant k. But 
it strongly depends on e being responsible for the energy 
that must be paid to tilt or roll two respective base pairs. 
Since a change of twist for constant Ri is proportional to 
a change in bond length the bond energy contributes to 
the twist persistence length explaining the dependence of 
Itw on k (compare Fig. EJ. 
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FIG. 8: (Color online) Comparison of analytical expressions 
Eqs. © and lllll for Ip and In (solid lines) with numerically 
calculated orientational correlation functions (data points) for 
2b = SA, fc = 64fcs^/A^ and e = 20, 60 [ksT] (from 
bottom to top). 



We also measured the mean-square end-to-end dis- 
tance (i?!;) and find that deviates from the usual 
WLC chain result due to the compressibility of the chain. 
So as to investigate the origin of the compressibility we 
calculate {R%) for the following geometry. We consider 
two base-pairs without spontaneous bending angles such 
that the end-to-end vector Re can be expressed as 



-Rb =^c, =^(Rit, -f Shb, -f Sin,). (29) 

i i 



The coordinate system {t^, b^, rii} is illustrated in Fig. ^ 
Ci denotes the center-center distance of two neighboring 
base-pairs. Since successive base-pair step parameters 
are independent of each other, and Ri and Sh and SI are 
uncorrelated the mean-square end-to-end distance {R%) 
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FIG. 10: (Color online) Comparison of the simulation data 
with e = 20kBT, k = 64fcsT/A^ 2b = llA, and T = 1, 2, 3, 5 
(from top to bottom) to Eqs. ll^ and ll^ (solid 

hnes). Using the measured bending persistence lengths and 
the stretching moduli we find a good agreement with the pre- 
dicted behavior. For T = 1 we obtain 7 — 6.O2A 

We compared the data for different temperatures T to 
Eq. l|5n|) using the measured bending persistence lengths 
Ip and stretching moduh 7 (see Fig. I10|l . The agree- 
ment is excellent. This indicates that transverse slide and 
shift fluctuations contribute to the longitudinal stretching 
modulus of the chain. 



FIG. 9: (Color online) Dependency of (a) bending persistence 
length Ip and (b) torsional persistence length l„ on the spring 
constant k, the width of the ellipsoids b and the energy depth 
e. We measured the persistence lengths for varying width sizes 
2b = 8, 9, 10, llA (red, green, blue, purple) and for two dif- 
ferent spring constants fc = 32 (plus), 64 (circles) [fesT/A ]. 
The bending persistence length depends solely on b and e. It 
gets larger for larger e and b values. But it does not depend 
on k (the curves for different k values corresponding to the 
same width b lie one upon the other). The torsional persis- 
tence length in turn depends on k, since a change of twist for 
constant Ri is proportional to a change in bond length. 



is given by 



«.2iV(m)^,-2/^(l-exp(-^ 
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(30) 



N denotes the number of base-pairs. Note that (SI) and 
(Sh) vanish. Using (cf) = (Ri2) + (Sh2) + (Sl^) the stretch- 
ing modulus 7 is simply given by 



(Ri; 



((Ri2)-(Ri)2) + (Sh2) + (Sl2) 



(31) 



C. Stretching 

Extension experiments on double-stranded B-DNA 
have shown that the overstretching transition occurs 
when the molecule is subjected to stretching forces of 
65pN or more . The DNA molecule thereby increases 
in length by a factor of 1.8 times the normal contour 
length. This overstretched DNA conformation is called 
S-DNA. The structure of S-DNA is still under discus- 
sion. First evidence of possible S-DNA conformations 
were provided by Lavery et al. 0, 0, [31 using atom- 
istic computer simulations. 

In principle one can imagine two possible scenarios how 
the transition from B-DNA to S-DNA occurs within our 
model. Either the chain untwists and unstacks resulting 
in an untwisted ladder with approximately 1.8 times the 
equilibrium length, or the chain untwists and the base- 
pairs slide against each other resulting in a skewed ladder 
with the same S-DNA length. The second scenario should 
be energetically favorable since it provides a possibility 
to partially conserve the stacking of successive base-pairs. 
In fact molecular modeling of the DNA stretching pro- 
cess 0, 0, 01 yielded both a conformation with strong 
inclination of base-pairs and an unwound ribbon depend- 
ing on which strand one pulls. 

We expect that the critical force fcrit where the struc- 
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FIG. 11: (Color online) Force-extension relation calculated 
by minimum energy calculation (black) and obtained by MC 
simulation (red) for 50 (red) and 500 (blue) base-pairs. The 
red solid line represents the analytical result of the WLC. (in- 
set) The deviation between energy minimization (black dotted 
line) and MC in the critical force is due to entropic contribu- 
tions. 



tural transition from B-DNA to overstretched S-DNA oc- 
curs depends only on the GB energy depth e controUing 
the stacking energy. So as a first step to find an appropri- 
ate value of e as input parameter for the MC simulation 
we minimize the Hamiltonian with an additional stretch- 
ing energy Epuii = fcis+i, where the stretching force 
acts along the center-of-mass axis, with respect to Ri, 
SI and Tw for a given pulling force /. Fig. 1111 shows 
the resulting stress-strain curve. First the pulling force 
acts solely against the stacking energy up to the crit- 
ical force where a jump from L{fcrit-)/ Lq ~ 1.05 to 

L{fcrit+)/Lo = VrF+W/Ri w 1.8 occurs, followed 
by another slow increase of the length caused by over- 
stretching the bonds. Lq = L{F = 0) = Ri denotes the 
stress- free center-of-mass distance. As already mentioned 
three local minima are obtained: (i) a stacked, twisted 
conformation, (ii) a skewed ladder, and (iii) an unwound 
helix. The strength of the applied stretching force de- 
termines which of the local minima becomes the global 
one. The global minimum for small stretching forces is 
determined to be the stacked, twisted conformation and 
the global minima for stretching forces larger than fcrit 
is found to be the skewed ladder. Therefore the broad- 
ness of the force plateau depends solely on the ratio of 
Z/Ri determined by the geometry of the base pairs S and 
the bond length I = 6.0A. A linear relationship is ob- 
tained between the critical force and the stacking energy 
e so that one can extrapolate to smaller e values to ex- 
tract the e value that reproduces the experimental value 
of fcrit ~ 65pN. This suggests a value of e « 7. 

The simulation results of the previous sections show 
several problems when this value of e is chosen. First of 
all it cannot produce the correct persistence lengths, the 
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FIG. 12: (Color online) (a) Probability distribution function 
of the center-center distance of successive base-pairs for / = 
(red), 140 (green), 200 (blue) pN. (b) Mean squared values 
of rise (red, plus), shift (green, crosses), slide (blue, stars), 
and center-of-mass distance (purple, squares) for neighboring 
base-pairs as a function of the stretching force /. The dashed 
line corresponds to the S-DNA center-of-mass distance. (Tw) 
of the resulting S-DNA conformation vanishes as predicted by 
Eq. ©. 

chain is far to flexible. Secondly the undistorted ground 
state is not a B-DNA anymore. The thermal fluctuations 
suffice to unstack and untwist the chain locally. That is 
why one has to choose larger e values even though the 
critical force is going to be overestimated. 

Therefore we choose the following way to fix the pa- 
rameter set {b,e,k}. First of all we choose a value for 
the stacking energy that reproduces correctly the persis- 
tence length. Afterwards the torsional persistence length 
is fixed to the experimentally known values by choos- 
ing an appropriate spring constant k. The depth of the 
base-pairs has also an influence on the persistence lengths 
of the chain. If the depth b is decreased larger fluctua- 
tions for all three rotational parameters are gained such 
that the persistence lengths get smaller. Furthermore 
the geometric structure and the behavior under pulling is 
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very sensitive to b. Too small values provoke non-B-DNA 
conformations or unphysical S-DNA conformations. We 
choose for b a value of llA for those reasons. For e = 20 
and A: = 64 a bending stiffness of Ip — 170bp and a tor- 
sional stiffness of Z„ — 270bp are obtained close to the 
experimental values. We use this parameter set to simu- 
late the corresponding stress-strain relation. 



The simulated stress-strain curves for 50 base-pairs 
show three different regimes (see Fig. Ill|) . (i) For small 
stretching forces the WLC behavior of the DNA in ad- 
dition with linear stretching elasticity of the backbones 
is recovered. This regime is completely determined by 
the chain length N. Due to the coarse-graining pro- 
cedure that provides analytic expressions of the persis- 
tence lengths depending on the base-pair parameters (see 
eqs. (|Hll,C3) it is not necessary to simulate a chain of 
a few thousand base-pairs. The stress-strain relation 
of the entropic and WLC stretching regime (small rel- 
ative extensions L/Lq and small forces) is known analyt- 
ically |20ll53| . Since we have parameterized the model in 
such a way that we recover the elastic properties of DNA 
on large length scales the simulation data for very long 
chains will follow the analytical result for small stretching 
forces, (ii) Around the critical force fcrit ~ 140pN which 
is mainly determined by the stacking energy of the base- 
pairs the structural transition from B-DNA to S-DNA 
occurs, (iii) For larger forces the bonds become over- 
stretched. Our MC simulations suggest a critical force 
fcrit ~ 140pN which is slightly smaller than the value 
fcrit ~ 180pN calculated by minimizing the energy. This 
is due to entropic contributions. 



In order to further characterize the B-to-S-transition 
we measured the mean values of rise, slide, shift, etc. as 
a function of the applied forces. The evaluation of the 
MC data shows that the mean values of shift, roll and 
tilt are completely independent of the applied stretching 
force and vanish for all /. Rise increases at the criti- 
cal force from the undisturbed value of 3.3A to approxi- 
mately 4.0A and decays subsequently to the undisturbed 
value. Quite interestingly the mean value of slide jumps 
from its undisturbed value of to ±5A (no direction is 
favored) and the twist changes at the critical force from 
7r/10 to 0. The calculation of the distribution function 
of the center-center distance c of two neighboring base- 
pairs for / = 140pN yields a double-peaked distribution 
(see Fig. I12|l indicating that part of the chain is in the 
B-form and part of the chain in the S-form. The contri- 
bution of the three translational degrees of freedom to the 
center-center distance c is shown in Fig. ^1 The S-DNA 
conformation is characterized by Ri = 3.3A, SI = ±5A 
and Tw — 0. In agreement with Refs. we ob- 

tain a conformation with highly inclined base-pairs still 
allowing for partial stacking of successive base-pairs. 
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FIG. 13: (Color online) Contour plot of rise [A] versus slide 
[A] for the S-DNA conformation. 

V. DISCUSSION 

We have introduced a simple model Hamiltonian de- 
scribing double-stranded DNA on the base-pair level. 
Due to the simplification of the force-field and, in par- 
ticular, the possibility of non-local MC moves our model 
provides access to much larger length scales than atom- 
istic simulations. For example, Ah on a AMD Athlon 
MP 2000-1- processor are sufhcient in order to generate 
1000 independent conformations for chains consisting of 
iV = 100 base-pairs. 

In the data analysis, the main emphasis was on de- 
riving the elastic constants on the elastic rod level from 
the analysis of thermal fluctuations of base-pair step pa- 
rameters. Assuming a twisted ladder as ground state 
conformation one can provide an analytical relationship 
between the persistence lengths and the local elastic con- 
stants given by eqs. (O, lfTT|l . Future work has 
to show, if it is possible to obtain suitable parameters 
for our mesoscopic model from a corresponding analysis 
of atomistic simulations |54| or quantum-chemical cal- 
culations js^. In the present paper, we have chosen a 
top-down approach, i.e. we try to reproduce the exper- 
imentally measured behavior of DNA on length scales 
beyond the base diameter. The analysis of the persis- 
tence lengths, the mean and mean squared values of all 
six base-pair parameters and the critical force, where the 
structural transition from B-DNA to S-DNA takes place, 
as a function of the model parameters {b, k, e} and the 
applied stretching force / suggests the following param- 
eter set: 

2b = llA (32) 
e = 20kBT (33) 

k = 64fcBr/A^ (34) 

It reproduces the correct persistence lengths for B- 
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DNA and entails the correct mean values of the base- 
pair step parameters known by X-ray diffraction studies. 
While the present model does not include the distinc- 
tion between the minor and major groove and suppresses 
all internal degrees of freedom of the base-pairs such as 
propellor twist, it nevertheless reproduces some experi- 
mentally observed features on the base-pair level. For 
example, the anisotropy of the bending angles (rolling is 
easier than tilting) is just a consequence of the plate-like 
shape of the base-pairs and the twist-stretch coupling is 
the result of the preferred stacking of neighboring base- 
pairs and the rigid backbones. 

The measured critical force is overestimated by a fac- 
tor of 2 and cannot be improved further by fine-tuning of 
the three free model parameters {b,k,e}. fcrit depends 
solely on the stacking energy value e that cannot be re- 
duced further. Otherwise neither the correct equilibrium 
structure of B-DNA nor the correct persistence lengths 
would be reproduced. Our model suggests a structure 
for S-DNA with highly inclined base-pairs so as to en- 
able at least partial base-pair stacking. This is in good 
agreement with results of atomistic B-DNA simulations 
by Lavery et al. 0, They found a force plateau 

of 140pN for freely rotating ends The mapping to 

the SOP model yields the following twist-stretch (Ri-Tw) 
coupling constant km^Tw = iC~^)m,Tw = 267/A. km^Tw 
is the microscopic coupling of rise and twist describing 
the untwisting of the chain due to an increase of rise 
(compare also Fig. EJ. 

Possible applications of the present model include the 
investigation of (i) the charge renormalization of the 
WLC elastic constants '56|, (ii) the microscopic origins 
of the cooperativity of the B-to-S transition and 
(iii) the influence of nicks in the sugar-phosphate back- 
bone on force-elongation curves. In particular, our model 
provides a physically sensible framework to study the in- 
tercalation of certain drugs or of ethidium bromide be- 
tween base pairs. The latter is a hydrophobic molecule 
of roughly the same size as the base-pairs that fluoresces 
green and likes to slip between two base-pairs forming 
an DNA-ethidium-bromide complex. The fluorescence 
properties allow to measure the persistence lengths of 
DNA ^. It was also used to argue that the force plateau 
is the result of a DNA conformational transition jll| . 

In the future, we plan to generalize our approach to a 
description on the base level which includes the possibil- 
ity of hydrogen-bond breaking between complementary 
bases along the lines of Ref. |30L l3ll |. A suitably pa- 
rameterized model allows a more detailed investigation 
of DNA unzipping experiments 58] as well as a direct 
comparison between the two mechanism currently dis- 
cussed for the B-to-S transition: the formation of skewed 
ladder conformations (as in the present paper) versus lo- 
cal denaturation [EH, [63, • Clearly, it is possible to 
study sequence-effects and even more refined models of 
DNA. For example, it is possible to mimic minor and 
major groove by bringing the backbones closer to one 
side of the ellipsoids without observing non-B-DNA like 



ground states. The relaxation of the internal degrees 
of freedom of the base-pairs characterized by another 
set of parameters (propeller twist, stagger, etc.) should 
help to reduce artifacts which are due to the ellipsoidal 
shape of the base-pairs. Sequence effects enter via the 
strength of the hydrogen bonds {Eqc = S.OfcsT versus 
Eat = l.SfcsT) as well as via base dependent stack- 
ing interactions 35]. For example, one finds for guanine 
a concentration of negative charge on the major-groove 
edge whereas for cytosine one finds a concentration of 
positive charge on the major-groove edge. For adanine 
and thymine instead there is no strong joint concentra- 
tion of partial charges . It is known that in a solution 
of water and ethanol where the hydrophobic effect is less 
dominant these partial charges cause GG/CC steps to 
adopt A- or C-forms 62^ by a negative slide and posi- 
tive roll motion and a positive slide motion respectively. 
Thus by varying the ratio of the strengths of the stack- 
ing versus the electrostatic energy it should be possible to 
study the transition from B-DNA to A-DNA and C-DNA 
respectively. 



VI. SUMMARY 

Inspired by the results of El Hassan and Calladine [s^l 
and of Hunter et al. "s?, 'sEj we have put forward the idea 
of constructing simplified DNA models on the base (-pair) 
level where discotic ellipsoids (whose stacking interac- 
tions are modeled via coarse-grained potentials js^H^) 
are linked to each other in such a way as to preserve the 
DNA geometry, its major mechanical degrees of freedom 
and the physical driving forces for the structure forma- 
tion [l|. 

In the present paper, we have used energy minimiza- 
tion and Monte Carlo simulations to study a simple repre- 
sentative of this class of DNA models with non-separable 
base-pairs. For a suitable choice of parameters we ob- 
tained a B-DNA like ground state as well as realistic 
values for the bend and twist persistence lengths. The 
latter were obtained by analyzing the thermal fiuctua- 
tions of long filaments as well as by a systematic coarse- 
graining from the stack-of-plates to the elastic rod level. 
In studying the response of DNA to external forces or 
torques, models of the present type are not restricted 
to the regime of small local deformations. Rather by 
specifying a physically motivated Hamiltonian for arbi- 
trary base- (step) parameters, our ansatz allows for real- 
istic local structural transitions. For the simple case of a 
stretching force we observed a transition from a twisted 
helix to a skewed ladder conformation. While our results 
suggest a similar structure for S-DNA as atomistic simu- 
lations JT], the DNA model studied in this paper can, of 
course, not be used to rule out the alternate possibility 
of local strand separations [s^ |63, |^ . 

In our opinion, the base(-pair) level provides a sen- 
sible compromise between conceptual simplicity, com- 
putational cost and degree of reality. Besides provid- 
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ing access to much larger scales than atomistic simula- 
tions, the derivation of such models from more micro- 
scopic considerations provides considerable insight. At 
the same time, they may serve to validate and unify an- 
alytical approaches aiming at (averaged) properties on 
larger scales [H IH EO, El E3 ■ Finally we note that the 
applicability of linked-ellipsoid models is not restricted to 
the base-pair level of DNA as the same techniques can, 
for example, also be used to study chromatin |63 . 163. l65| . 
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